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THE KLEIN-GORDON EQUATION ON Z^ AND THE QUANTUM 

HARMONIC LATTICE 

VITA BOROVYK AND MICHAEL GOLDBERG 

Abstract. The discrete Klein-Gordon equation on a two-dimensional square lattice satisfies an 

£^ I— ^ £°° dispersive bound with polynomial decay rate Itl"^'**. We determine the shape of the light 

(.^ ■ cone for any choice of the mass parameter and relative propagation speeds along the two coordinate 

#vj ' axes. Fundamental solutions experience the least dispersion along four caustic lines interior to the 

, light cone rather than along its boundary, and decay exponentially to arbitrary order outside the 

Oj' light cone. The overall geometry of the propagation pattern and its associated dispersive bounds 

are independent of the particular choice of parameters. In particular there is no bifurcation of the 

number or type of caustics that are present. 
CO ' The discrete Klein-Gordon equation is a classical analogue of the quantum harmonic lattice. In 

Cn I the quantum setting, commutators of time-shifted observables experience the same decay rates as 

the corresponding Klein-Gordon solutions, which depend in turn on the relative location of the 
^ ' observables' support sets. 

< 

c^ ; 1. Introduction 

The wave equation uu — An = on M?'^^ is explicitly solved via Poisson's formula, in which 
initial data u{x,0) = g{x), ut{x,0) = h{x) determines the unique solution 

>' , sign(t) /■ h{y) + \{g{y)+Vg{y)-{y-x)) ^ 

U{X,t) — ^ j r-7y j TTy "2/ 

2vr J\y-x\<\t\ ^t^-\y-x\^ 



m 
^^ . 

^' • at any t ^ 0. More generally the Klein-Gordon equation uu — An -|- m?u = with the same initial 

data has the solution 



^. u{x,t) = ^-^ / lh{y) + -ig{y) +Vg{y) ■ {y - x)))^^^=JL=^ dy. 

When 771 = the two equations coincide. It is clear in these formulas that the propagator kernel 
. . , is radially symmetric, and that all information from the initial data travels at finite speed. Both 

r^ • equations also possess a well-known dispersive property that \u{x, t)\ < C|t|~^' ^ provided the initial 

C^ . data are sufficiently smooth and decay at infinity. 

This paper investigates the propagation patterns and dispersive bounds for solutions to a family 
of discrete Klein-Gordon equations on Z'^ xM}, 

uu{x,t) — X],=iAj(n(x -|- ej,t) + u{x — ej,t) — 2u{x,t)^ + uj'^u{x,t) = 

Hu 

n(x,0) = g{x) 
ut{x,0) = h{x) 

with w, Ai, A2 > fixed parameters. The conserved energy is given by 

, 2 

(1.2) E{t) = - ^ {ul{x,t)+uj'^u^{x,t) + ^\j{u{x + ej,t)-u{x,t)f) 

x£7? j = l 



(1.1) 
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and (by analogy with the continuous setting) the values of \j suggest propagation speeds of a/Aj 
along their respective coordinate directions. We have chosen the letter lo for the mass parameter 
in order to highlight connections between (jl.ip and the quantum harmonic lattice system. 

The one-dimensional wave equation provides a certain degree of inspiration; when A = 1 its 
fundamental solution is expressed in terms of the Bessel functions J\y_x\{t) (see |llj). There are 
three main asymptotic regimes. For |i| ^> |y — x| there is oscillation with amplitude |i|~^'^. When 
1^1 ^ |y ~ 2;| the propagator is nonzero (hence there is some rapid transfer of information) with 
exponential decay at spatial infinity on the order of ((2|y — x|)~^et)l^~^L For |y — x| = \t\ + 0{t^'^) 
the propagator kernel reaches its maximum size of approximately |t|~^'^. This bound is most easily 
obtained by applying van der Corput's lemma to the Fourier representation of J\y_x\{t) . 

Unfortunately the discrete wave and Klein-Gordon equations in higher dimensions do not sepa- 
rate variables as does the Schrodinger equation. The fundamental solution can still be determined 
as a superposition of plane waves, with size bounds in the different regimes resulting from station- 
ary phase principles. The two and three-dimensional isotropic wave equations (Aj = 1, a; = 0) were 
analyzed by Schultz [16], with a curious set of outcomes. In addition to the expected wavefront 
expanding radially at ly — x| ~ |t|, there is a secondary region of reduced dispersion traveling at 
somewhat lower speed. In two dimensions the region lies along an astroid-shaped curve with diam- 
eter \/2 l^h i^ three dimensions the region follows a cusped and pointed surface of a similar nature. 
Surprisingly, some global dispersive bounds are dominated by behavior when |y — x| belongs to the 
secondary set, even though this occurs well inside the overall propagation pattern. 

In the two-dimensional discrete Klein-Gordon equation, each plane wave Uk{x) := e satisfies 
Huk = 7^(A:)ufc, with 7^(fc) = w^ + ^ ■ 2Aj(l — cos kj) and k ranging over the fundamental domain 
[— 7r,7r]^. The solution of (jl.ip is given formally by 



u{x,t) = cos{tvH)g -\ -j= — h and ut{x,t) = vH sm{tvH)g + cos{tvH)h. 



The operators involved act on a plane wave Uk by multiplication by cos{t^{k)), ^^"^21 , and 
7(/c) sin(t7(A;)), hence the fundamental solutions of (II. ip in physical space will be the inverse Fourier 
transform of those three functions. For all practical purposes these are oscillatory integrals over 
the torus k S [— vr, vr]^ of the form given in (12. 2p . whose asymptotic behavior is governed by critical 
points of the phase function t^{k) ziz x ■ k. 

Three distinct regimes again emerge: critical points are absent for x ^ t and exponential decay 
is observed by following the analytic continuation of 7(A;) into {[— '7r,7r] -|- zR}^. Quantitative ex- 
ponential bounds are given in Theorems 12.61 and 12.71 For generic values of {x,t) inside the "light 
cone" (i.e. x = tV^{k) for some k), stationary phase arguments lead to a bound of |t|~^. Along 
the boundary of the light cone, and within the secondary region introduced above, degenerate 
stationary phase estimates yield polynomial time decay with a fractionally smaller exponent. For 
fixed i 7^ there is a global bound of order |t|~^'^ with maxima occuring near the four cusps of the 
astroid curve. This is a faster rate of decay than the discrete Schrodinger equation on 7?, where 
separation of variables leads to a |i|~^'^ bound instead. Further details about the structure of the 
Klein-Gordon propagators are summarized in Theorem 12.31 and its corollaries. 

The problem of generalizing van der Corput's lemma to two and higher dimensions has a long 
history in harmonic analysis. Since the degeneracy that arises during our computations is not 
very severe, we are able to follow Varchenko's 1976 exposition |17j . Modern techniques for the 
general resolution of singularities may be needed for applications where the domain has a more 
intricate periodic structure than Z^, and especially in high-dimensional settings. In those cases one 
may employ methods and results by Greenblatt [6] in two dimensions or Collins, Greenleaf, and 
Pramanik [4] in higher dimensions. At one point in Corollary 13.41 we also invoke a recent result by 
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Ikromov and Miiller [7] regarding the stability of degenerate integrals under linear perturbation of 
the phase. 

Unlike in the continuous setting, varying parameters Ai, A2 is not equivalent to performing a 
diagonal linear transformation on x because the domain Tl? and its Fourier dual both lack a dilation 
symmetry. We show in this paper that the dispersion pattern for (jl.ip retains the same topological 
and geometric structure found in [16] for all values of a;,Aj. In particular there is no choice of 
parameters that generates exceptional degeneracy of the dispersive estimate or bifurcation of the 
dynamical system. Separately we show that interactions outside the light cone are all subject to 
an exponential bound, and that exponential bounds of any desired order are achieved by setting 
\y — x\l\t\ sufficiently large. The latter statement is akin to (and readily implies) Lieb- Robinson 
bounds (cf. [8], [H], [3], [5], [l5], [l2], [H]) for the corresponding quantum harmonic lattice. The 
background and details of this application are presented in Section [3 

Section [2] enumerates the precise statements and illustrations of our main results. The proof is 
sketched out in the following section, assuming a number of propositions about the critical points 
of t7(A;) — x-k. Section [4] takes a detailed look at the Taylor series expansions of 7 in order to verify 
those assertions. The concluding section translates our results about the discrete Klein-Gordon 
equation into dynamical properties of the quantum harmonic lattice. 

2. Main results 

In this section we state the result describing the long-time behavior of the solution of (jl.ip . We 
start with estimates of a slightly more general oscillatory integral and our main result is based on 
these estimates. 

For the function 

(2.1) 7(A;)= a;2 + ^2Aj-(l-cosA;j-) , /e = (fci, ^2) G [-vr, ^]2, 

introduce 

(2.2) I{t, X, 11) = -^ / e^C^— *^W)r?(A;) dk, 

where x € 1? , t € M, and r/ is a smooth test function on [— vr, vr]^ with periodic boundary conditions. 
The asymptotic behavior of oscillatory integrals is generally influenced by local considerations, in 
which case r/(x) may be assumed to have compact support in a fundamental domain of 'M? /2tt1? 
(for example as part of a partition of unity). In that case ry can be extended by zero to a function 
on all of M^ and one may define 

(2.3) /(t, X, ,?) = -!-/■ e^C^ -*7W)^(A;) dfc, 

{2ttY Jf,2 

for all x G M^ up to a unimodular constant whose value is exactly 1 if x E Z^. It is often convenient 
to consider x of the form x = vt, with v a fixed vector in R^ (representing velocity), so that the 
integral ()2.3p can be written as 

(2.4) i(t,vt,ri) = -^ f e'^^-^''^ri(k)dk, where <pJk) := k ■ v - j(k). 

[2iry J^2 

Estimates on I{t, x, rj) will follow from restricting the corresponding bound on (j2.4p to examples 
with X = vt &T?. 

The long-time behavior of (j2.4p is dictated by the highest level of degeneracy of the phase within 
the support of r/. In the absence of critical points, integration by parts multiple times yields a rapid 
(faster than polynomial) decay. In the case where critical points are present but are non-degenerate, 
the standard stationary phase argument provides \t\~^ decay for the integral. Finally, if there are 



4 VITA BOROVYK AND MICHAEL GOLDBERG 

degenerate critical points, the decay is slower and more careful analysis is needed to determine its 
exact order. It is easy to see that for any point k* G [— vr, tt]^, there is a choice of the velocity v such 
that k* is a critical point of (pvik) (namely, v = 'Vj{k*)). The order of degeneracy of the phase at 
that point is determined by second and higher-order derivatives of 7, as the linear component is 
cancelled by subtracting k* -v. We introduce a partition of [— vr, vr]^ with respect to the degeneracy 
order of 7, 



(2.5) 



[-TT, nf = Ki\jK2LI K3 



where 



(2.6) 



Ki = {k e [-vr,7r] 
K2 = {ke [-vr,^]2 

K3 = {ke [-TT.TTf 



2 -.det D^-f{k) /O}, 
■.detD'^{k) = 0, (e.V)S(A:)/0}, 
: det D'^-f{k) = 0, (e • V)S(^) = 0}. 



In the definition of K2 and K3, C stands for an eigenvector of the 2x2 matrix D'^'y{k) corresponding 
to the zero eigenvalue. 

The analysis of Section H] allows us to describe the structure of this partition in detail. Propo- 
sition UT] notes in particular that the rank of D^^{k) is never zero, so the direction of ^ is always 
well defined. 



1,2,3, defined in (12. 6p possess the following properties. 



Lemma 2.1. The sets Ki, ■ 

• Ki consists of four points related by mirror symmetry across the coordinate axes. 

• K2 consists of two closed curves, one around the origin and the other around the point 
(7r,7r), with the four points of K^ removed from the latter curve. 

• Ki = [— 7r,7r]^\ {K3UK2)- This set consists of three open regions: the interior of the small 
closed curve around zero, the interior of the closed curve around the point (vr,7r), and the 
area of the compactified torus enclosed between these two curves. 

The structure of the partition is displayed in FigureU\ (similar to Figure\S^. 



-vr 



IT 



& 



TT 



-^ ki 



-TT 



Figure 1. Sets Ki, K2, K3 



Proof. The statements follow immediately from equation ()4.5p and Corollary 14.8 



D 



THE KLEIN-GORDON EQUATION ON T? AND THE QUANTUM HARMONIC LATTICE 5 

To see the connection between the velocities and possible degeneracies of i;^„, let us analyze 
images of the sets Ki in the velocity space: 

V3 = {v G M : there exists k G i^3 such that v = V7(A;)}, 
(2.7) V2 = {v £R'^\V3: there exists k £ K2 such that v = Vj{k)}, 

Vi = {v £R^\ {V2 U V3) : there exists k £ Ki such that v = V-yik)}, 
Vo = {veR^ : for all k e [-7r,7r]2, v / V7(A:)}. 

-7r,7r]2^]R2 defined by 
V{k)=Vj{k), 



Alternatively, under the mapping V 

(2.8) 

sets (j2.7p admit the representation 



(2.9) 



V3 = V(i^3), 
V2 = V{K2), 

Vi=ViKi)\{V3UV2) 
yo = M2\V([-^,7r]2). 



The collection {V^}f=o forms a partition of R?. More detailed description of its structure is provided 
in the following result. 

Proposition 2.2. Let the sets {Vi}f^Q be defined by (j2.7p . Then they are located as shown on 
Figure [H There are two simple closed continuous curves ^1 and ^2 around the origin that split 
the plane into three open regions. More precisely, ^1 encloses a convex region and ^2 consists of 
four concave arcs that meet at cusps. The four vertices of these cusps form the set V3. The union 
of ^1 and ^2, with the four points that belong to V3 removed, is V2. The union of the two inner 
bounded open regions is Vi. The outer-most unbounded region is Vq which has boundary ^1. 




^.Xl 



Figure 2. Sets Vb, ¥l, V2, and V3. 



We can now state the main result describing the connection between a velocity v and the decay 
order of an oscillatory integral (|2.4p with phase function 7. 
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Theorem 2.3. Let the sets {Vi]^^^ be defined by (j2.7p . t] be a smooth periodic function on [—it, tt]^, 
and integral I{t,x,r]) defined by (|2.2p . Then for any fixed 6 > 0, there exist constants Cq,Ci, C2 
and C3 depending on r] such that 

C3 



(2.10) for all x with dist (x, tV^) < t6, we have \I{t, x,ri)\ < 



t|3/4' 

C2 



(2.11) for all x with dist (xjtVs) > t6 and dist {x,tV2) < t6, we have \I{t,x,i])\ < 



t|5/6' 

c 

(2.12) for all x with dist {x,t{Vs U V2)) > tS and dist(2;,iyi) < t6, we have \I{t,x,r])\ < -r-j, 

C 

(2.13) given N > 1, for all x with dist(x,t(V3 U V2 U Vi)) > 5, we have \I{t,x,r])\ < y-n^- 

Each of Co, Ci, and C2 depend on 6 (and Cq also depends on N). The values of Cq - C2 are 
expected to approach infinity as 5 approaches zero, while C3 is independent of 6. 

If one is only concerned with the worst possible decay among ah x G Z^ , the simpler statement 
is as follows. 

Corollary 2.4. Let the integral L{t,x,ri) be defined by (|2.2|) . Then there exists C = C{rj) > such 
that for all x G Z^, 

(2.14) |/(t,x,7?)|<^. 

Corollary 2.5. As a special case of Theorem \2.'A the propagators of the discrete Klein-Gordon 
equation (jl.lj) are recovered by choosing rj = 7™, m = —1,0,1. Hence the values of u{x,t) and 
ut{x,t) for a solution of (|l.ip both satisfy (J2.10p - (j2.13p . provided the initial data g, h are supported 
at the origin. The result extends to all g, h supported in B{0,R) by superposition, once enough 
time has elapsed that 5t > 2R. 

The values of the exponents in (j2.10p - ()2.12p are dictated by the worst degeneracy degree of 
critical points of the phase function. For example, if a velocity belongs to the set Vi and is 
relatively far from V2 and V3, all the critical points of (p^ will be uniformly non-degenerate. In this 
case the decay rate of an oscillatory integral is |t|~'^''^ in arbitrary dimension d (hence it is |t|~^ in 
dimension two). 

Let us now briefly describe how the rates produced by velocities that are near V2 and V3 are 
computed (see Section \3A\ for details). If u G V2 U V3, then (j)^ has at least one degenerate critical 
point k* G [— TT, vr]^. Then the Taylor series expansion of (p^ near its critical point takes the form 

(2.15) (f,,(k) = k ■ V - ^(k) = CO + Y. Cn,m{kl-KT{k2-kir. 

n,m>0 
n+m>2 

This Taylor series is said to be supported on the set of indices (m, n) where Cm,n 7^ 0. Roughly 
speaking, one computes the leading-order decay of the corresponding oscillatory integral by mea- 
suring the distance from the origin to the convex hull of the Taylor series support, then taking the 
reciprocal. However the support is not invariant under changes of coordinates, so one must first 
choose an "adapted" coordinate system that maximizes this distance |17j . It turns out that for 
each function cj)v{k) the linear coordinate system that diagonalizes the Hessian matrix is adapted 
(Lemma 14.91 verifies this property in the one case where it is not readily apparent). According to 
the definitions ()2.6p and ()2.7p . the Taylor series of (pv with respect to these coordinates has more 
vanishing low-order terms if w G V3 as compared to f G V2. Therefore its Newton polyhedron 
lies further away from the origin, and the oscillatory integral decays more slowly. The Newton 
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polyhedra associated with v £ V2 and f G V3 are sketched in Figures [3] and H] and give rise to the 
exponents in (|2.1ip and (|2.10p respectively. 

If < dist(f,V2 U V3) < 6, then (j)^ does not have a degenerate critical point itself, but it is 
related to the degenerate phase functions described above by a small linear perturbation. The fact 
that oscillatory integral estimates are stable under such perturbations is proved in [7]. 

According to its definition, Vq C M^ consists of velocities that produce phase functions with no 
critical points in the domain of integration. As a result one can recover polynomial decay (in x) 
of I{t, X, 1]) of any order by repeated integrations by parts. In fact for solutions of the discrete 
Klein-Gordon equation the decay satisfies a number of exponential bounds. 

Theorem 2.6. For every /i > there exists constants < v^ < -(1 + 2\/Ai + A2 sinh(/x/2)) and 
C^ < uj + 2v'Ai + A2 cosh(;u/2) such that 



(2.16) 
(2.17) 
(2.18) 



1 



(27r)2 
1 



(27r)2 



(27r)2 



cos{t-f{k))e^^-'' dk 
sin(t7(fc)) ^^fc.^. ^^ 

-l{k)sin{t^{k))e'^-''dk 



< e 



< e 



-^fi{\x\~v^\t\) 



-fj-{\x\-'"M) 



< C e~'^^'''''~''f''*'^ 



The upper bound for u^ as stated in Theorem 12.61 behaves as expected for large fx (see Corol- 
lary 2.2 in [12]) but it has some evident drawbacks over the rest of the range. First, the sharp value 
of Vfj^ must be an increasing function of fx so the apparent asymptote as ;U — )• is an artifact of the 
calculation. In addition the estimates (j2.16p and (j2.17p don't show any time-decay when applied 
to a point x G tVo with -^^ < f^. The last result shows that in fact every x G tVo is subject to an 
effective exponential bound. 

Theorem 2.7. Let Vq be the set defined in (j2.7p . For any x £ M? with j £ Vq there exists fi > 
and constants Ci < 00, C2 < \/u}'^ + 4(Ai -|- A2) such that 

1 



(2.19) 
(2.20) 
(2.21) 



(27r)2 
1 



cos{tj{k))e'''-''dk 



<e 



- Iidist{x , tVi) 



{2ny 



sin(t7(A:)) ^.fc.^ 
.12 j{k) 



(2vr)2 



^{k)sm{t-f{k))e'^-''dk 



< (^^g-Mdist(x,tVi) 



3. Proof of the main results 



3.1. Proof of Theorem 12.31 The material of this section is presented in the following order: we 
start with some background inforination on oscillatory integrals, followed by the main local results 
(Lemma 13.11 and Corollary 13. 4p , and the proof of Theorem 12.31 is then obtained from local estimates 
through a partition of unity argument. 

Consider an oscillatory integral in several variables 



(3.1) 



Iit,v) 



jt,pik) 



7]{k)dk, 



where r] is supported in a neighborhood of an isolated critical point k* of (j) (we follow the notation 
introduced in J17j and [6]). When it is convenient to do so, one may apply an affine translation so 
that the critical point is located at the origin. If supp(r7) is small enough, 7(t, rj) has an asymptotic 
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expansion 

oo 

(3.2) /(i,r?)?«e**'^(°)^(dj(r/)+4(r7)ln(t))t-"j, as t ^ oo, 

i=o 

where Sj is an increasing arithmetic progression of positive rational numbers independent of i]. The 
oscillatory index of the function (j) at k* is defined to be the leading-order exponent sq. We assume 
that So is chosen to be minimal such that in any sufficiently small neighborhood U containing k* 
either do{rj) or dgC^) i^ nonzero for some r/ supported in U. 

Estimates ()2.10p and (12. lip are essentially statements about the oscillatory index of (t)v{k) at its 
critical points for different values of v. The following algorithm assists in their computation. 

Suppose (j) is analytic with a critical point at k* . Locally there is a Taylor series expansion 
4>{k) = Co + ^ Cn{k — A;*)", with the sum ranging over all n G Z^ with ni + . . . + n^ > 2. Let 
K C Z"*" be the collection of all indices n for which c„, 7^ 0. 

Newton's polyhedron associated to (p at its critical point k* is defined as the convex hull of the 
set 

|J(n + M^), 

neK 

where W^ is the positive octant {x G M'^ : Xj > for I < j < d}. We denote Newton's polyhedron 
of (p by N^{(p). Newton's diagram of (j) is the union of all compact faces of N^{(j)). Finally, the 
Newton distance d{4>) is defined as (i(^) = inf{i : {t,t) G A^+(0)}. 

Note that the vanishing of Taylor coefficients is affected by changes to the underlying coordinates, 
thus each local coordinate system y generates its own sets K^ and N^{(p) and Newton distance 
dy{(j)). Define the height of an analytic function (p at its critical point to be h{(f>) := sup{d^((/>)}, 
with the supremum taken over all local coordinate systems y. A coordinate system y is called 
adapted if #(<^) = h{(p). 

It was shown in [17] (p. 177, Theorem 0.6) that under some natural assumptions, the oscillation 
index sq of a function 4' is equal to l/h{4>). 

We now compute the height of the phase of the integral I(t, rr, rj) at its critical point(s). Recall 
that with the notation v = j, the phase of I{t,x,r]) can be written in the form (I2.15p . Given a 
point k* £ [— 7r,7r]^, the value v = \7j{k*) € M? is the unique choice for which cpy has a critical 
point at k* . Denote the height of this (f)^ at k* by h{k*). 

Lemma 3.1. The height function h{k*), defined above, is constant on each of the sets {Ki}^^^ 
defined in (j2.6p with the following values: 

(1) ifk* G Ki, thenh{k*) = 1, 

(2) ifk* £ K2, then h{k*) = 6/5, 

(3) ifk* G A'3, then h{k*) = 4/3. 

Proof. Note that ^{k) and <pv{k) differ by a linear function, so their derivatives coincide except 
at the first order. In the simpler case k* G i^i, detD^7(A;*) = dei D"^ (^^{k*) 7^ by definition. 
Moreover, the determinant of the Hessian of (pv at k* remains non-zero in any local coordinate 
system, thus h{k*) = d/2 = 1. 

Next consider k* G K2. In this case detD^7(fc*) = and (^ • V)'^7(A:*) 7^ 0, where ^ is an 
eigenvector of D^"f{k*) corresponding to the zero eigenvalue. The mixed second-order derivative 
vanishes because D^"f{k*) has orthogonal eigenvectors, and by Proposition 14. II it is guaranteed that 
(^ ■ V)^7(A;*) 7^ 0. This information suffices to compute the Newton's distance of ^^ at k* in the 
linear coordinate system with axes {C , C} ^'Hd given by coordinates k — k* = y\^ + y2^- 

The associated Newton's polyhedron is of the form displayed on Figure El with the Newton's 
distance being 6/5. 

In order to show that {2/1,2/2} is an adapted coordinate system, and therefore h{k*) = 6/5, we 
use the following result from Varchenko: 
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n2 . 




Figure 3. Newton's polyhedron and Newton's distance of the Taylor series corre- 
sponding to ^-u for k* G K2 

Proposition 3.2 ( |171 part 2 of Proposition 0.7]). Assume that for a given series f = Yl^nU"', the 
point {d{f),d{f)) lies on a closed compact face T of the Newton's polyhedron. Let aiui + a2r^2 = "m 
be the equation of the straight line on which T lies, where ai, 02, and m are integers and ai and 02 
are relatively prime. Then the coordinate system y is adapted if both numbers ai and 02 are larger 
than 1. 

The Newton polyhedron displayed on Figure [3] has only one compact face (that also contains the 
point {d{<j)y) , d{(j)y))) , which lies on the line with the equation 3ni + 2n2 = 6. Since 2 and 3 are 
relatively prime, the coordinate system is adapted by Proposition 13.21 

Finally, let k G K^. To provide a more concise notation, let d^ and dt± indicate the directional 
derivatives ^ • V and ^-^ ■ V respectively. These also serve as partial derivatives dy^ and dy-^ with 
respect to coordinates {yi,y2}- By definition of K3 we have 

(3.3) dl^/ik*) = 0, d^d^±j{k*) = 0, 9|x7(^"*) / 0, and dlj{k*) = 0. 

On the other hand, one can show that df'^{k*) and d'ldp^'^{k*) do not both vanish (see Lemma [4. 9p . 
The possible Newton's polyhedra that arise are indicated in Figure [H however, the Newton distance 
is equal to 4/3 in all situations: 

n2 . 



a. 



-N{(t>,)- 




ni 



Figure 4. Possible Newton's polyhedra and Newton's distance of the Taylor series 
corresponding to (p^, for k* G K^ 

Moreover, in all situations the face containing the point {d{(j)y) , d{(j)y)) lies on the line with the 
equation 2ni + n2 = 4, and Proposition 13.21 does not apply. To verify that the system {yi,y2} is 
adapted in this case, we use a different result from Varchenko: 

Proposition 3.3 ( [n\ Proposition 0.8]). Assume that for a given series f = Ylcny"', the point 
{d{f),d{f)) lies on a closed compact face T of the Newton's polyhedron. Let airii + n2 = m be the 
equation of the straight line on which T lies, where ai and m are integers. Let 

(3.4) /p(y) = ^c„y" and P{yi) = fr{yi,l). 

ngr 
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If the polynomial P does not have a real root of multiplicity larger than ?n(l + 02) ^ , then y is a 
coordinate system adapted to f. 



For the face F, displayed on Figure HI we have 



(3.5) 



My) 



Pivi 



dl^k* 



OlMk* 



2 ^ 5|5^x7(fc* 



-yt + 



2 dhik* 

■yiyi+ ^ 



-yl + ^-^ yi + 



24 

dh{k*) 



yt 



'5 



24 



The discriminant of P is 
(3.6) V 



a|a^x7(F)y dl^ik*) dl^j{k*] 



24 



1 
12 



3{dld^.-f{k*)f - dl^^{k*)dl^{k* 



It follows from Lemma [4.9l that this discriminant is nonzero whenever k* G K^,, thus P can have real 
roots of multiplicity at most one. On the other hand, m(l + 02)"^ = 4/3, and by Proposition 13.31 
the coordinate system {yi,y2} is adapted. D 

As was stated earlier, the height of the phase function determines the decay order of an oscillatory 
integral in a neighborhood of its critical point. In [T7] it is shown that the oscillation index of a 
phase 4) is equal to l/h{(p), giving both upper and lower bounds for the decay rate. More recently, 
Ikromov and Miiller in [^ showed that the upper bound is stable under linear perturbations of the 
phase function. Their result, combined with Lemma l3.H brings us the following 

Corollary 3.4. Let sets {-fCjjf^x be defined by ()2.6p and fix k* € [— vr, vr]^. Then there exist a 
neighborhood of k* , Vtk*, o^nd a positive constant Ck* such that for all rj supported in ilk*, 



(3.7) 

(3.8) 

(3.9) 

for all X € 



i(t,x,7j) 
I{t,x,r]) 
i{t,x,r]) 



< Cfc. \\v 

< Ck" ||r/||c3 

< Ck* Ihlica 



1 



'■) 



|i|3/4' 

1 
^576' 
1 

w 



tfk*GK3, 

ifk*eK2, 

ifk* £ Ki. 



Proof. By a direct consequence of Theorem 1.1 in [7], for a point k* G 
neighborhood Vtk* and a positive constant Ck* such that 



-7r,7rl , there exist a 



(3.10) 

for all X G 
claim. 



5**(<^"W+^-'=)r?(A;)dA; 



fe* 



<Ci 



'=* ll^ll U|h(fc*)' 



P and r/ supported in il^*. This, together with the result of Lemma |3. 11 proves the 

D 



We can extract some additional important information about the neighborhoods Vt^* , introduced 
in Corollary 13.41 that will be useful for the main part of the proof. Specifically, if k* G Ki then 
the corresponding neighborhood Ofc. does not contain any points from K2 U K^, and if k* G K2 
then Ofc* is disjoint from K^,. To prove the latter claim, suppose there is a point /cq S K^ that also 
belongs to O^*. Then with v = V7(/i:o) the oscillation index of (f)^ in a neighborhood of feo is equal 
to 3/4. For ry supported in a small neighborhood of /cq inside of ilfc*, and x = tv, the asymptotic 
lower bound dictated by (13. 2p and statement (3) of Lemma l3 . 1 1 contradicts the decay rates of (13. 8p 
and (13:91) . 

At last, we need the following well-known estimate. 
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Lemma 3.5. //suppr/ does not contain any critical points of the phase function k ■ x — t'y{k) then 
for any M > 0, 

(3.11) |/(t,x,7?)|<C(M,r/,d). ^ 



\t\M ■ 

where d is the infimum of \x/t — V7(A;)| over the support of 7]. 

Proof of Theorem \2.3[ Fix S > 0. The nonstationary phase bound (j2.13p follows immediately from 
the construction, as the gradient of x ■ k — t'y{k) must have magnitude at least dist(x,tVi). 

To prove (j2.10p . Take the system of neighborhoods {Qk}ke[-n.TT]^ described in Corollary 13. 4[ By 
construction {^k}ke\-n,n]'^ covers [— 7r,7r]^ and we can choose a finite sub-cover, say {^j}j^i- Now, 
let a collection of smooth functions {uoj},^-^ form a partition of unity with respect to {$7^} •^-^, then 

(3.12) \Iit,x,rj)\ <Y,\I{t,x,7ij)\ =Y,\Ht,x,7ij)\, 

j J 

where rjj = rjUj, is supported in Oj. Since for t away from zero every integral that satisfies (j3.8p 
or (j3.9p also satisfies (|3.7p . and all three are uniformly bounded for all times, we have 

(3-13) \nt,x,r^)\ <Y,Cj ||r/,||c;3(M2) ^ = C{v)i^. 

Note that even though (j3.13p holds for all j; G M^, better estimates are available when x is removed 
from tVs. 

To prove the estimates (j2.1ip and ()2.12p we need to refine our construction of the cover so that 
the |t|~^''^ bound in (j3.7p is never invoked (in the latter case one should also avoid applying (|3.8p ). 
The following construction will suit both situations. 

The function V, defined by (j2.8p . is uniformly continuous on [— 7r,7r]^, so we can choose an 
< e = e{6) < 7r/2 such that 

diam(V(Se)) < 6/2 
for every ball B^ of radius e. At each k G [— 7r,7r]^ define a smaller (^-dependent neighborhood 

(3.14) nkiS) = nknB,{k), 

where O^ is again as in Corollarv 13.41 As before, pick a finite sub-collection of {^fc('^)}fcgf-7r,7ri2 
that is also a cover of [— vr, vr]^, say {^kj{6)}jLi, and generate a partition of unity loj subordinate 
to this cover. For simplicity of notation we will write Qj = Qk (6) with j £ {1, 2, . . . , N}. 

Sort the neighborhoods Qj according to the location of their "center" point kj. For each ni = 
1, 2, 3 let Jm '■= {j € {1, ... , N} : kj £ K^}, where K^n are the sets defined in (j2.6p . The discussion 
following Corollary 13.41 indicates that 

(3.15) ( IJ % j n Kg = and ( [J % j n A's = 0. 

Suppose X G Z^ is chosen so that dist(x,tV3) > tb. In other words, |x/t — V7(A;*)| > b for any 
k* G K-^. Moreover |x/t — V7(A;)| > 5/2 for all k G U/eJs ^J because each neighborhood has radius 
at most e. Split the sum (|3.12p into two parts 

(3.16) |/(*,x,r?)| < Y, |/(t,x,r/,-)| + Y. l^^'^'^jOI- 

iG-^3 j<tJ-A 

Lemma 13.51 applies to each term in the first sum, with d = 6/2. Terms in the second sum are 
bounded by (j3.8|) or p.9p . The slowest time-decay out of these has the rate |t|~^'^ from (13. Sp . 
which verifies (|2.1ip . 
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The argument is similar if dist(x,i(V3 U V2)) > t^- One splits (j3.12p in the parts 
(3.17) 



\I{t,x,rj)\< Y. \i{t,x,r,j)\+ Yl \^(t,x,i]j)\, 



and once again Lemma 13.51 applies to each term in the first sum, with d = 5/2, and terms in 
the second sum are bounded by (j3.9p . This is sufficient to verify ()2.12p . completing the proof of 
Theorem [231 D 

3.2. Proof of Exponential Bounds. 

Proof of Theorem \2.6[ Note that 7^(fc) extends to a complex-analytic function on k £ C'^ that is 
periodic under the shifts kj — )• kj + 27r, j = 1,2. After composition with the holomorphic map 
cos(ty^), the same is true of cos(t7(A:)). By shifting the contour of integration for ki and k2, the 
left-hand quantity in (j2.16p is equal to 



-fj.\x\ 



(2vr)^ 



cos{tj{k + ifiA))e"'''' dk 



< max 

fcG[-7r,7r]2 

< max e 



cos{t'y{k + ifi-^))\e 



W 



-fj.\x\ 



Imtyik+ifj.T-r) 



-MFl 



where v^ = fi~^ max{|Im7(/c -|- ik)\ : k £ [— vr, vr]^, |A;| = fi}. Referring back to the definition of j{k) 
in (|2.ip . one obtains a bound v^ < j^V^i + -^2 sinh(/_f/2) by applying the inequality |Im y/z'^ + w'^\ < 
y^(Imz)^ -|- (Imu;)^ for pairs of complex numbers. 

The same argument applies to the sine propagator as well, thanks to the bound ^^ < g|im2|_ 
By shifting the integration contour as above, the left-hand quantity in (j2.17p is equal to 



-f^\x\ 



sin{tJ{k + ^^lJ^)) 



(27r)2 7[_^^^]2 t-/{k + in^A 



dk 



< t max e 



\lrat'y{k+ifi 



w 



"MF 



< e 



-^J.(\x\~v^,\t\) 



for any v^> 11 ^(1 -|- max{|Im7(A; -|- ik)\ : k G [— 7r,7r]^, \k\ = /i}). 

The computation for (|2.18p is essentially identical to (j2.16p except that it contains an extra factor 
of max{|7(A; -|- ik) : \k\ = fi}, estimated here by w -|- 2-v/Ai -|- A2 cosh(/i/2). 

D 



Proof of Theorem 2.1 . It will be important to note that Vq is the complement of a convex subset 
of the plane. This is stated as part of Proposition 12.21 and will be proved in Section [H 

For each /2 G R^, shifting contours of integration into the complex plane leads to the bound 



(2vr)^ 



cos(t7(A:))e*^'dA; 



-pt-x 



(2vr) 



cos{tj{k + iil))e"'''' dk 



< max el^'^*^(^"+^^)le-'^-^ 



By assumption j lies outside the convex balanced compact set {V7(A;) : k G [— tt, vr]^} = Vi. 
The complex derivative of 7 indicates that Im7(A: -|- ijl) = {V^{k)) ■ jl + o(|/i|), and the implicit 
constant in o(|/i|) converges uniformly across k G [— 7r,7r]^. Choose fi > small enough so that 



|lm7(A: + ijl) - fl ■ {Vj{k))\ < -dist(f , ViM 



THE KLEIN-GORDON EQUATION ON I? AND THE QUANTUM HARMONIC LATTICE 13 

whenever |/i| = 2^1 . Then 



(27r 



^ ' cos(t7(A;))e*^(iA; 



■12 



< mf max f.\{i^^m■-A^^^^^{^^y^)i^f.-^■^^ 

/i|=2/^ fce[— 7r,7r]2 

= mf max e(*V7(fc)-)-Aedist(x,tyOM 

|/i|=2/^ fce[-7r,7r]2 
— 2/xdist(x,iVi) dist(a::,tyi)^ 

The first equahty follows from the fact that ^^{k) is an odd function, so the absolute value can be 
optimized with either sign. The second equality asserts a geometric principle that given a closed 
convex set S and a point x ^ S, 

sup {y — x) ■ V > —\v\ dist(x, S) 
yes 

with equality taking place ii y £ S minimizes the distance and v is parallel to x — y. The argument 
for the sine propagator is essentially identical as the extra factor of t that also appears in the proof 
of Theorem 12.61 can be overcome by choosing a slightly smaller value of /i > and introducing 
a large constant Ci. The value of C2 is limited by estimating the maximum of \j{k + ik)\ over 

4. Properties of the phase function 
Throughout this section we will be using the variables (a, b) defined by 

(4.1) a{k) = cos ki, b{k) =cosk2, 

along with k = {ki,k2)- The function 7, defined in (|2.ip . as a function of (a, 6) takes the form 

(4.2) 7(a,6) = (w2 + 2Al(l-a) + 2A2(l-fe))^^^ 

A straightforward calculation shows that the Hessian matrix of 7 can be written in the form 

. . 2 n\- ^ /^ Aia7^(fc) - Af (1 - g^) -A1A2 sin /ci sinfe2 

[^.6) u ^[K) - ^3^^^ (^ -AiA2sinA;isinA:2 \2hl'ik) - Ai(l - \?) 

and, thus, 

detZ52^(A:) = 4^ [ah-i^ih) - Ai6(l - a^) - A2a(l - \?)) 

^ [k) ' 

(4.4) = -^ [ahu? - AiKl - af - A2a(l - hf) . 

Proposition 4.1. For every choice ofuj, Ai,A2 > 0, there is no k £ [— 7r,7r]^ such that D'^^{k) is 
the zero matrix. 

Proof. The off-diagonal entries vanish only if sin fci = or sin /c2 = 0. Without loss of generality, 
suppose sin/ci = 0. Then a = coski = ±1, so that 9| ^{k) = ±Ai7~"'^(A;) 7^ 0. If sinA:2 = 0, then 
9| 7 (A;) is nonzero for similar reasons. D 

One of our goals is to describe the set of zeros of the Hessian determinant, 

(4.5) $1 = {A; G [-vr, irf : det Z?2^(A:) = 0}. 

However, it will be convenient to first study the zeros of det D"^^ as a function of (a, 6): 

(4.6) ri = {(a,6) G [-1,1]^ : det 1)27(0, 6) =0}. 
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Using the notation 

F{a, b) = abj^{a, h) - Ai6(l - a^) - A2a(l - 6^) 

(4.7) = ahu"^ - Ai6(l - af - \2a{l - bf, 

we have that det D'^^{a, 6) = if and only if F(a, b) = (since Ai, A2, 7(A;) 7^ 0), and therefore, 

(4.8) ri = {(a,6)G[-l,l]2:F(a,6) = 0}. 

Sometimes it will be convenient to treat i^ as a function of k, and in those cases we will keep the 
same notation, F = F{k). Note also that with the notation (14. 7p the Hessian matrix takes the 
form: 

1 / XidbF — A1A2 sin/ci sin/c2 

X2daF 



(4.9) 



D^-f{k) 



73(fc) 



-A1A2 sin/ci sin /c2 



Lemma 4.2. The equation F{a,b) = defines an implicit function b = Bp{a) in [—1, 1]^ that is 
locally continuously differentiable at any (a, 6) G Fi, |6| ^ 1. It has the following properties: 

(1) For any a G [—1, 1], there exists at most one b G [—1, 1] so that F{a, b) = 0. 

(2) For any (a, 6) Gri\{(0,0)}, |6| / 1, 



(4.10) 



(4.11) 



dBp 
da 



ia,b) 



Ai 6^(l-a^) 
'A2a2(l-62) 



<0, 



and 



^C'")-^' 



(3) The set Fi defined in (j4.8p is the graph ofb = Bp{a), which consists of two continuous arcs 
T\ U F^ as displayed in Figure\M The first arc, T\, is located in the first quadrant and is 
convex. The second arc, T^, passes through the second and fourth quadrant and is concave. 

The equation F{a,b) = also defines an implicit function a = Apib). All of the statements in 
this lemma remain true if (a, Ap, Ai) and {b,BF,X2) switch roles. 





b 

1 






-1 


— ^ 


V 


1 






\ 






-1 







Figure 5. Set Fi 



Proof. For fixed a G [—1,1] the function F{a,b) is quadratic with respect to b and has at most 
one solution in the interval 5 G [—1, 1]. The same is true if b is held fixed and one seeks the value 
a = AF{b) for which F{AF{b),b) = 0. As a result Bp is a well-defined function over some subset 
of [—1, 1] and Ap serves as its inverse. 
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One can write out the value of Bp{a) explicitly using the quadratic formula and derive all the 
stated properties from this expression. We present a more general approach here in preparation for 
subsequent computations where an exact formula is not readily available. 

The function F is continuously differentiable and 



db 



^— (a, b) = au^ - Ai(l - af + 2A2a(l - b) = - (abuj^ - Ai6(l - af) + 2A2a(l - b) 

1 



^ ,F(a, b) + A2a(l - 6)^) + 2A2a(l - b). 



Therefore, 



(4.12) ^^("' ^) = ^2a(l - b^) for all (a, b) G Ti 



and this quantity is nonzero so long as a 7^ and \b\ < 1. When a = one can compute directly 
dF 
db 



that If (0,6) = -Ai / 0. Therefore 



dF 



= if and only if |6| = 1. 



An identical argument applied to the variable a shows that 

r)F 

(4.13) a^-(a, b) = XMl - a^) for all (a, b) G Ti 

oa 

and furthermore that ^(o;6)|p = if and only if \a\ = 1. 

In order to prove equation (j4.10p . we differentiate F{a, 6) = implicitly with respect to a. Taking 
advantage of ()4.12p and ()4.13p the resulting expression reduces to 



db b/adaF{a,b)\ Xib'^{l-a- 



^ • ^ da a\bdbF{a,b)J A2a2(l-62) 

for all (a, b) G Fi away from the origin. At the origin, ()4.1ip is obtained directly from the facts 
that ^(a, 0) = — A2 and ^(0,6) = — Ai. This is consistent with the implicit derivative in ()4.10p 
since both statements demand that the ratio 6^/a^ converges to (A2/Ai)^ as (a, 6) approaches the 
origin along Ti. 

By definition Fi must be the graph of Bp, which is continuously differentiable with negative 
slope whereves it lies inside (—1, 1)^ and has slope zero when \a\ = 1. Given that < Bp(—1) < 
Bf{1) < 1 and < Ap{—1) < Ap{l) < 1, it follows that Fi consists of two separate arcs. One 
arc, denoted by T\, connects the points {Ap{l),l) and {l,Bp{l)) within the first quadrant. The 
second arc, denoted by Ff, connects (— 1, i?i;'(— 1)) to {Ap{—1),—1) and passes through the origin 
(since -F(0, 0) = 0) along the way. 

Finally, a routine derivation shows that 



S 



0\ /i2 



which is clearly positive if a, 6 > 0, thus proving that F} is convex. Using again that F(a, 6) = 0, 
we rewrite the second derivative in the form 

2A 6^ 

= -A^ a4(l-62)3 "^ (^' + ^i(^ - ")'(2 + ") + ^2(1 - b)\2 + 6)) 
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So long as < |a|, |6| < 1, the sign of this second derivative is determined by the sign of ab, which 
is negative everywhere on Tf except the origin. A separate calculation shows that 



(4.16) 



d^ 



(0,0) 



-2A2(u;^ + 2Ai + 2A2) 
A? 



<0, 



finishing the proof. The above expression can be simplified further by noting that (x! +2A1+2A2 = 7 
when a = b = 0. D 

Remark 4.3. Recalling that a = cos/ci and b = cos ^2 for points {ki,k2) G [— vTjTt]^, we can recon- 
struct the graph of $1 (defined in ()4.5p ) from the graph of Fi (see Figure[6]). The arc T\ corresponds 
to the closed curve around zero. The origin in the a6-plane has the four points {±ir/2, ±7r/2) as its 
inverse image and the arc Vf turns into the closed curve around the point (tt, tt) on the compactified 
torus [— 7r,7rp. 



-vr 



^2 i 



TT 



-TT 



& 



TT 



-* h 



Figure 6. Set <l>i 

By definition a function (j)v{k) = k ■ v — ^{k) may have a degenerate critical point only at 
k* € <I>i, and this occurs with the choice v = \/^{k*). For each k G <l>i, let ^ = ^{k) be an 
eigenvector of D'^^{k) corresponding to its zero eigenvalue (this is unique up to scalar mutiplication 
by Proposition 14. ip . It follows that d't^{k) = 0, where the notation d^ = ^-V indicates a directional 
derivative as in Section [3j According to the partition (j2.6p . points A: G $1 belong to either K2 or 
K^ depending on whether the third derivative of 7 in the direction of ^ is also zero. The following 
result is helpful. 

Lemma 4.4. Let [/ C M"' &e a neighborhood of a point k^ and let f G C^{U). Assume that the 
Hessian matrix D^f has a zero eigenvalue of multiplicity one at ko and let (^ be a corresponding 
eigenvector. Then 

(4.17) difiko) = 
if and only if 

(4.18) d^{detD^f){ko) = 0. 

Proof. Apply a unitary change of variables to change the coordinate system to one that diagonalizes 
the matrix (D^/)(A;o) and in which ^ points in the direction of ei. In the new system, the only 
non-zero term of the gradient of det(L'^/)(/c) at ko is the gradient of (9ii/)(fe) at ko multiplied by 
a nonzero scalar - the product of all nonzero eigenvalues of {D'^f){ko). On the other hand, we have 
that 

1 ^T „0 „ /, N J, 1 



(4.19) 



idtjm 



m 



,eD'f{k)i 



m 



2 5|/(^)> 



THE KLEIN-GORDON EQUATION ON I? AND THE QUANTUM HARMONIC LATTICE 17 

showing that (V(det D^/))(/i:o) is a non-zero scalar multiple of V(9?/)(A;o). D 

Lemma 14.41 allows us to determine whether points A; G $1 satisfy df^{k) = by identifying the 
set of solutions of 

(4.20) d^{detD^-f){k) = 0. 
Using equations ()4.4p . ()4.7p . and notation (j4.ip we have 

(4.21) VdetD%(A;)|,^^^ = ^VFik) = -^(daFsinki, dbFsink^). 



The components of ^ can be constructed from the elements of the matrix (j4.9p . with one possible 
choice being 

Note that ^(/c) vanishes as sinfci approaches zero along $1 (for example by applying (j4.13p ). and 
in fact 

(4.23) -]—^(k)^ f, ° 1 /O 

smki \X1smk2J 

when ki — t- 0, itvr along this curve. The combination of ()4.22p with ()4.2ip yields 

(4.24) a^(detD%)(A;) = -^i^g^((5.F)2 + Ai(l-62)a,F), k G <^i. 

One should not be concerned with the vanishing of sin ki in this formula as it can be counteracted by 
modifying ()4.22p by a suitable scalar multiple. Vanishing of the second factor determines whether 
A; e $1 belongs to K2 or K3. Using (j4.12p and ()4.13p . we have 

aH {{daFf + Ai(l - b')dbF) = Xlb'il - a^f + \1\2a\l - b'f 

and thus if A; G <I>i (equivalently if (a, b) G Fi), then ()4.20p holds only if 

(4.25) G{a, h) = Ai63(l _ a^f + \^a\l - b^f = 0. 

The function G is symmetric in a and b and is rather elegant, but it turns out not to be ideal for 
our purposes. Restricting our view to A; G $1, we are again free to add any multiple of F to G and 
work with that object instead. We therefore introduce 

(4.26) G{a,b) = G{a,b) + aH^F{a,b) 

= u^a^b^ + Aife^(l - 3a^ + 2a^) + A2a^(l - 36^ + 26^). 

The function G maintains the property that among all A; G <^i, ()4.20p holds only if G{a, b) = 0. We 
will describe some features of the set 

(4.27) Fa = {(a, b) G [-1, 1]^ : G{a, b) = 0} 

as an independent object before seeking out its intersection with Fi. The following result is an 
analog of Lemma 14.21 and we provide it for the sake of completeness. 

Lemma 4.5. The set F2 defined in (j4.27p is nonempty and is of the form displayed in Figure^ it 
consists of two continuous arcs, F2 = F2UF2. The equation G{a,b) = defines an implicit function 
in [—1,1]^ that is locally continuously differentiable with respect to a at any (a, 6) G F2, |6| 7^ 1. 
In particular, the arc F2 represents the graph of a function which we denote by b = Bcia). The 
following properties hold: 
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(1) For any {a,b) G ri\{(0,0)}, |6| / 1, 



(4.28) 

and 

(4.29) 



dBg 
da 



ia,b) 



Ai b\l - g^) 
■A^a4(l-62) 



<0, 



da 



(0,0) 



1/3 



(2) The first arc, T^, passes through the third quadrant. The second arc, r^, is located in the 
second and fourth quadrant. 

The arc F^ is also the graph of a function a = Acib). All of the statements (or their analogs) in 
this lemma remain true if {a,AG,Xi) and {b,BG,X2) switch roles. 





b 








1 








\ 






-1 


\ 




1 






^- 






-1 







Figure 7. Set r2 



Proof. The proof is largely identical to that of Lemma 14.21 and we omit the common details. One 

f (0,0) and f 



difference is that a slope at the origin cannot be determined from the ratio of ^(0, 0) and ^(0, 0) 



as both quantities are zero already. 

Note that G(0, 5) = Ai6^ is zero only if 6 = 0. When a 7^ 0, write r(o) = ^ to obtain the 
expression 



(4.30) 



G{a, r) = a3((Air3 + A2) + a\-3Xir^ - SAsr^) + a^{Lo^ + 2Ai + 2A2)r3). 



For a fixed value of a, the solutions of G{a, r) = occur at the roots of a cubic polynomial whose 
coefficients depend smoothly on a. When a = the polynomial is Xir^ + A2, which has a single 
transversal root at r = — (^)^'^. The implicit function theorem provides a neighborhood of a = 
and a continuous function r{a) along which G{a, r{a)) = 0. 

By definition (0, 0) = lim r(a) = — (t^)"*^'"^- Once again the result is consistent with the gen- 

da a->o "^1 

eral implicit derivative (|4.28p because the ratio (6^/a^) converges to (A2/Ai)^/^ as (a, 6) approaches 

the origin along F2. 

D 

We are most interested in the intersection points of Fi and F2, which are described in the 
following result. 



Lemma 4.6. Let the curves Fi and F2 be defined by ()4.8p and (14.270 . respectively. Then 
(i) i/Ai < A2, FinF2 = {(0,0), (a*, 6*)}, with a* < < b* , 
(ii) i/Ai > A2, FinF2 = {(0,0),(a*,r)}, with b* <0<a*, 



(4.32) !li^ (a) = !1^ (a)^, for all (a, 6) G Fi n Ta, (a, 6) / (0, 0). 
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(iii) i/Ai = A2, rinr2 = {(o,o)}. 

In the last case, we will use the notation {a*,b*) = (0,0). 

Proof. Consider the case Ai < A2- The proof in the case Ai > A2 is identical. 

Note that since the origin belongs to both Fi and r2, it is obviously in their intersection. 
Lemma g^l and Lemma |13] show that L^ n Li = L} n r2 = and thus, Li n r2 = L^ n Lf . 
Next, according to Lemma 14.21 ^1 is concave, and with the assumption A2 > Ai its slope at zero is 
less than —1 (see equation (j4.1ip ). As a result, for all {a,b) G Tf in the fourth quadrant, |6| > \a\, 
except for the origin. Define 

(4.31) a+ =max{a : (a,6) Grinr2}. 

If we assume that Tf intersects Fg in the fourth quadrant away from the origin, we have that 
< a~^ < a, where a = Ap{—1). Comparing formulas (|4.10p and ()4.28p gives 

dBc , , dBp , ,b^ 
da da a^ 

Both derivatives are negative, and furthermore |6| > \a\ along this part of the curve F^. Conse- 
quently 

and thus, Bcia) < BF{a) for all a G (a^,a^ + e) for some small e > 0. On the other hand, it follows 
from the results of Lemma 14.51 that Bcid) > — 1 = Bpid), implying that there must exist another 
point a' G [a"*" + e, o) such that Bcia') = Bpia'). This contradicts the definition of a+ and we may 
conclude that Fi n F2 = {(0,0)} in the fourth quadrant. 

We claim that Ff intersects F2 exactly once in the second quadrant away from the origin. First, 
since 

(^•^^) ^(°)<^(»)' 

there is a small e > so that Bp{a) > Bda) for all a G (— e, 0). However for the value a = ^^(1), 
one has Bp^d) < 1 = Bcid). The Intermediate Value Theorem implies the existence of a' G (a,0) 
such that Bp{a') = Bq^o'), giving rise to at least one non-origin point of intersection of F| and F| 
in the second quadrant. 

Note that G{a,-a) = -{oj'^ + AXi)a^ + {X2- Xi)a^{l-3a'^ -2a^) < over the interval a G [-1,0). 
Together with ()4.29p that implies that F2 lies above the line b = —a within the second quadrant, 
and it follows from (li:32|) that > ^(a) > ^(a) whenever Tj and F^ intersect with a < 0. 
On the other hand, if there were multiple points of intersection, the orientation of crossing would 
have alternating signs. One concludes that Ff n F2 contains a single point (a*, 6*) in the second 
quadrant along with the origin. 

Now consider the case Ai = A2 = A, and let (oq, 60) G Fi n F2. Then G(ao, &o) = 0, where G is 
defined in (|4.25p . Since Ai and A2 are equal, G admits the factorization 

(4.35) G{a, b) = X{a + b){a^ + b^ + ab{a^b^ - 2ab - 1)). 

The second factor is zero if and only if \a\ = \b\ = 1, however those points do not belong to Fi. 
Hence, 60 = —o-o- Plugging this into (|4.7p . one can see that when Ai = A2, F{ao, —ao) = if and 
only if ao = 0. The intersection of Fi and F2 at the origin is not transversal in this case, but instead 
the two curves are tangent without crossing. D 

Remark 4.7. Introduce the set K* as the pre-image of the point (a*, 6*) under the map (|4.ip . The 
set K* consists of four points on the set $1 that are located as shown on Figure [8j Note that in 
the case Ai = A2, K* = {(±f ,±f)}. 
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Figure 8. Set ^i along with the four-point set K* 



The following is an easy consequence of Lemma 14.61 

Corollary 4.8. Assume that det D'^-){k) = and let ^ = ^{k) be an eigenvector of D'^^{k) corre- 
sponding to the zero eigenvalue. Let K* be the set defined in Remark \4-T\ Then dc^{k) = if and 
only if k £ K* . 

Proof. First recall that the assumption det D^'y{k) = is equivalent to the fact that the corre- 
sponding (a, 6) G Fi (see (I4.6P ). 

Assume that d^^{k) = 0. In Lemma 14.41 we showed that under our main assumption this 
condition is equivalent to (j4.20p . This, in turn, implies that G{a,b) = (or (a, 6) G F2), where 
{a,b) is an image of k under the map (j4.ip . We may therefore conclude that {a,b) £ Fi n F2. 
According to Lemma [^^ in the case Ai = A2, {a,b) = {a*,b*) and thus, k S K* . If Ai 7^ A2, (0,6) 
is either (a*, 6*) or the origin. However, a direct inspection (see equation (|4.24p ) shows that when 
(a,6) = (0,0), 

(4.36) d^idet D^^){k) = ±^ (A? - A^) / 0, 



"7^W 



and we can exclude the origin from consideration. Thus, {a,b) = {a*,b*) and k E K* . 

The proof in the reverse direction is similar. Assume k £ K*, then it corresponds to (a*, 6*). 
Again, in the case Ai = A2, {a*,b*) = (0,0), and by K36\\ . d^ {det D'^-f){k) = a|7(fc) = 0. If 
-^1 7^-^2 5 both a* and b* are different from zero. Comparing equations ()4.24p - ()4.26p . we see that 
G{a*,b*) is a nonzero multiple of a* {b*)^d^ {det D^-f){k). Since G{a*,b*) = and a*{b*)^ / 0, we 
conclude d^{det D'^-f){k) = (9|7(A:) = 0. 

D 

Lemma 4.9. Let 7 be defined by (|4.2p and let k* G K* , where K* is defined in Remark \4-7\ 
Furthermore, let S, = ^{k*) be the eigenvector of the Hessian matrix ofj at k* corresponding to the 
zero eigenvalue and .^ be a vector orthogonal to ^ and of the same magnitude. Then 

(4.37) (a|7a|x7 - 3 [dp^^^f] {k*) + 0. 



<0. 



Proof. To prove equation (I4.37p . it is enough to show that 

(4.38) (5|7a|x7-2(5|a^x7)')(A;* 
The above expression admits the following short representation that we will use in our calculations, 

(4.39) (9|75|x7 - 2 (519^x7)') {k*) = UW' (9|detD%) (F). 
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Let us first prove (j4.39p . Indeed, we can re- write det 1)^7 in the new coordinates as 

(4.40) UfdetD^j = dl^dl^j - {d^d^±j)\ 
Differentiating tlie above equation we obtain 

At the point k*, the quantities 9|7, d^dc±^, and df^ all vanish, thus the second, third and fifth 
term of the above equation vanish as well, proving (|4.39p . 
Next, it is easy to see that 

(4.41) idldetD^,)ikn = {dlF)ikn^y 

where F is defined in (|4.7p . Finally, a direct calculation shows that in the case Ai 7^ A2, 

a*6*(l-(a*)2)(l-(6*)2) 



(4.42) (a|F) (F) = 2IICII V(^ 

and in the case Ai = A2 = A, 

(4.43) (a|F) (k*) = 
finishing the proof. 



(a*)2(l-(6*)2) + (6*)2(i_| 
-||ef(cu2 + 2(Ai + A2))<0, 



,*\2> 



<0, 



D 



4.1. Proof of Proposition 12.21 The curves of <l>i consist of points where det D'^'y{k) = 0, which 
are also the points where the "velocity map" V{k) = V7(A;) does not satisfy the hypotheses of the 
inverse function theorem. As a result the boundary of V([— 7r,7r]2) must be a subset of ^1 U ^2 = 
V{^i) as defined in Proposition 12. 2i 

Recall from Remark 14.31 that $1 has one closed curve around zero and a second closed curve 
around the point (vr,7r). Let ^1 be the image of the former under the velocity map and let ^2 be 
the image of the latter. The analysis of ^1 is more straightforward because the points of K* are 
not involved. 

In vector form, the velocity map is 

1 /Aisin/ci^ 
j{k) \A2sin/i;2^ 

Thus points /c in a given "quadrant" of the torus [— vr, vr]^ are mapped to the same quadrant of M^. 
The tangent line to <^i always points in the direction normal to V det D'^j{k), which by (j4.2ip 

is also normal to \/F[k). Suppose k travels along <^i with instantaneous velocity 



Vj{k) 



daF sinki^ 

Then k follows either loop of <I>i through the four quadrants of the compactified torus in order, and 
V{k) must wind once around the origin. 

At a local level, the differential dV{k) is the Hessian matrix D'^j{k), whose image is spanned by 
the direction ^^. Plugging (|4.9p into the Leibniz rule we determine that V{k) moves with velocity 

1 / XidhF — A1A2 sin/ci sinA;2\ f —d^F smk2^ 

-A1A2 sin/ci sinfc2 ^2daF J y daFsinki 

-XiidbF)^ - AiA2(l - a^)daFsmk2 



73 (A:) 



1 



(4.44) 



73 (A;) 
A1A2 

A1A2 

^3(fe) 



XiX2{l - b'^)dbF sin ki + X2{daFf sin k2 
Xib^il - a2)2 + A2a3(l - 62)2] f_a^-^,^^ 



G{a,b) 



-f sin/c2 
sin/ci 



b 
sinfci 
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Identities ()4.12p and ()4.13p are used multiple times between the second and third lines. 

The prefactor of AiA27~'^(fc) is positive for all k. The factor of G/{a^b) is strictly positive as 
k traces out the loop of $i circling the origin because a, 6 > here. The discussion leading up 
to Lemma 14.61 and Corollarv 14.81 shows that for general k G <l>i, the value of G/{o?b) changes sign 
when k crosses a point of K^ and at no other time. 

The vector with components (— | sin A:2, sin ki) points in the direction of ^^ and does not vanish 
while /c S $1 (the points where a = 6 = are handled by (j4.1ip ). Indeed one could choose this 
vector as the definition of ^ {k). Consider ^ as measured in polar coordinates. The path of V{k) 
turns to the left or the right depending on whether the polar angle of ^ is increasing or decreasing 
with k. The direction of change for this angle in turn depends on the sign of the determinant 

, /-fsinA;2 - (A2 "^K^ + Ai ^^) sin fci \ 1, , 2n2 . \ u^ 2n2\ 

det ". \ a\i-b^) ' = 12(^20(1-6) +Ai6(l-a) . 

y smfci -A2— ^-5 — ^smfc2 / " 

The left column of the 2x2 matrix above is ^^{k). The right column is its rate of change as k travels 

/ sin fci sin k^ _«\ / —dhF Smk'2\ 

along <l>i at the prescribed velocity, computed via the product I ^ n^ ] { Fi tp ' h 1 ' 

In the analysis of (j4.15p . this quantity is shown to have the same sign as ab, which is positive on 
the loop of $1 corresponding to 'I'l and negative on the loop corresponding to ^2- 

One last note is that the sign of ^ is constant on the connected components of ^i, so the direction 
of ^ winds exactly once around the origin on each loop. Combined with the preceding facts, it 
follows that 'I'l is a simple closed curve enclosing a convex region in M^, and that ^2 is a simple 
closed curve composed of four arcs with the opposite curvature from ^1. These arcs intersect at 
cusps located at the points V3 corresponding to values of A; G i^3 where G/{a?b) changes sign. 

It remains to be shown that ^1 is the boundary of V([— vr, vr]^) and that ^1, ^2 are disjoint. 
By comparing supplementary angles, it is clear that the extreme values of V{k) in any given 
direction must occur while |/ci|,|A;2| < ^. With the exception of k = (ib^,ib^), all points where 
V{k) G ^2 satisfy cos/ci cos A;2 < 0, so one of the coordinates is necessarily greater than ^. When 
1^1 1 = 1^2! = f 1 the vector ^^ which spans the image of DV happens to be collinear with V(A;), so 
these choices for k do not generate extreme points of V([— 7r,7r]^) in their respective directions. By 
process of elimination, the boundary must consist of ^1 alone. 

5. Applications to quantum systems 

In this section we apply the main integral estimates from Section [2] to obtain dispersive estimates 
in a class of infinite-volume harmonic systems. The material of the two introductory subsections 
follow [1], we provide it here for completeness. For more details, see [1] and references therein. 

5.1. Harmonic evolutions in finite volume. We first introduce a class of finite volume harmonic 
systems on two-dimensional square lattices. For an integer L > 1, denote A^ = {—L,L] C 1?. We 
associate with each x G A^ the position and momentum operators on L^(M, dq^): Qx is the operator 
of multiplication by Qx and px = —i^^- 

We denote by Qx and px the extensions of the corresponding single-site operators to the full 
Hilbert space 



(5.1) nL= (^ L'{R,dqx), 

xgA_l 

defined by setting 

(5.2) g^ = ]l(g)---(g)ll(g)gj,(g)]l(g)---(g)ll and p^^ = 1 (g) • • • (g) 1 (g) -i 



dqx 
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Operators qx and "Px are self-adjoint on Jii^ and satisfy the canonical commutation relations (CCR): 



(5-3) \Px,Vy\ 

for all x,y ^ Kl. 
The Hamiltonian 



(5.4) Hl 



and [qx.Py] 



■t:^ pI + w^ ql + y Xj {qx 



^^x^y^ : 



xGAt 



^ Qx + Z^ -J 



,)' 



is a self-adjoint operator on Hilbert space Hl and represents a system of coupled harmonic oscilla- 
tors. In the above, {ej}'^^^ represents the canonical basis vectors in Z^, the parameters Xj,uj > 
represent the coupling strength and the on-site energy. Finally, the Hamiltonian Hl is equipped 
with periodic boundary conditions, i.e., if a; G A^ but x + Cj A/,, set qx+Sj = Qx-(2L-i)e ■ 

Let B{Hl) denote the space of bounded linear operators over the Hilbert space T-Ll- We will refer 
to elements of B{T-Ll) as observables. As a self-adjoint operator, the Hamiltonian Hl generates the 
time evolution r/" on B{T-Ll), given by 



(5.5) 



T^{A) = e'^"^Ae~'*"^ for all t G M and aU A G B{nL) 



r/' is a well-defined, one parameter group of *-automorphisms, called the Heisenberg dynamics. 

A special class of observables called Weyl operators are important for the transition to the 
infinite-volume systems. For any / : A^ — )• C, a unitary operator 



(5.6) 



W{f) = exp 



Y, iMfix)]Qx + Mfix)]Px) 



XGAt 



is called a Weyl operator or a Weyl observable. One can show that 

(5.7) W{fr = W{-f) 
and 

(5.8) W{f)W{g) = e-^Im[</'9)]/2;y(/ + g) 

hold for aU f,g : Al ^ C. In addition, W{0) = 1 and \\W{f) - 1|| = 2 for each / / 0. An 
important fact is that the set of Weyl operators is invariant with respect to r/': there exists a 
family of operators T^^ : ^{Kl) -^ I^^{Al) such that 

(5.9) Tt{W{f)) = W{Ttf). 

Formula (j5.9p is verified by diagonalizing Hl with Fourier- type operators (cf. [I3]). The operators 
T/" have an explicit construction 



(5.10) 



tI^ = {U + V)F-^MtF{U* - V*) . 



— : X G Al \ and Mt is 



the operator of multiplication by e^*"^*, where 
(5.11 



l{k) 



\ 



OJ" 



+ A^Xj sirr^{kj/2) , /c G A^. 



Operators U and V , known as Bogoliubov transformations (cf. [9]), are given by 



(5.12) 



U = -F-^My+F and V = -F~^Mr_FJ, 
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where J is complex conjugation, and A/p^ is the operator of multipHcation by 

(5.13) r±(A;) = ^l=±v/^, 

V7(^) 

with 7 (A;) as in (j5.1ip . 

5.2. Harmonic evolutions in infinite volume. We start our discussion of infinite harmonic 
lattice with a review of general Weyl algebra formalism (see [10] and [2] for details). 

Let P denote a real linear space equipped with a symplectic, non-degenerate bilinear form a. The 
Weyl algebra over P, which we will denote by W(P), is then defined to be a C*-algebra generated 
by Weyl operators, i.e., non-zero elements W{f), associated to each f € T>, which satisfy 

(5.14) W{f)* = W{-f) for each / G P , 
and 

(5.15) W{f)W{g) = e-^'"(^'3)/2^(/ + g) for all f,g€V. 

It is well known (cf. [2], Theorem 5.2.8.) that such an algebra with additional properties that 
W{0) = 1, W{f) is unitary for all / G P, and \\W{f) - M = 2 for all f £V\ {0} is unique up to 
a *-isomorphism. 

In our case a convenient choice of V will be V = (?{^) or V = £^(Z^) with the symplectic form 

(5.16) a{f,g)=lm[{f,g)] iorf,geV. 

We denote the corresponding Weyl algebra by >V(^^(Z^)) or W(^^(Z^)). 

Let us return to the general case. Another result of Theorem 5.2.8 of j2] is that a group of real 
linear mappings {Tt}t<=M., Tt : V ^ V such that 

(5.17) a{Ttf, Ttg) = a{f, g) , for ah t G M, 

generates a unique one-parameter group of *-automorphisms r^ on W(P), such that 

(5.18) Tt{W{f)) = W{Ttf) foraU/GP. 

In order to use this result to introduce a harmonic dynamics on W(^^(Z^)), define Tj on l'^{l?) 
by (compare with (I5.10p ) 

(5.19) Tt = {U + V)F~^MtF{U* - V*) . 

Here F : £^(Z^) — )• L^([— vr, vr]^) is the unitary Fourier transform and Mf is the operator of mul- 
tiplication on L^([— 7r,7r]^) again by e^*'*'*, where the function 7 is defined as in (|2.ip . on [— 7r,7r]^ 
(compare to (j5.1ip ). Operators U and V are also defined as in ()5.12p with the appropriately 
extended objects. 

It is easy to see that {Tt}t£R is a family of real linear mappings, that satisfies the group properties 
To = 1, Ts+t = Tg oTt. Moreover, for each fixed t, Tt is symplectic, i.e., 

(5.20) lm[{Ttf,Ttg)]=lm[{f,g)], 

therefore (I5.17P is satisfied (see [T] for details) . Thus we may conclude that there exists of a unique 
one-parameter group of *- automorphisms on W(^^(Z^)), denoted by rj, such that 

(5.21) rt{W{f)) = W{Ttf) for aU / G £2(Z2) . 

The family Tt is called the infinite volume harmonic dynamics on W(^^(Z^)). 
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5.3. Main Results. We start with a standard estimate on the norm of the commutator of two 
Weyl observables. Using the Weyl relations (|5.15p . we get 



{Tt{WU)),W{g)\ = {W{Ttf)-W{g)W{Ttf)W{-g)}W{g) 



(5.22) 



Jlm[{Ttf,g)] 



} W{Ttf)W{g) . 



Since all the Weyl operators are unitary, we have 



(5.23) 



\\[rtiWif)),Wig)]\\ 



1 



JIm[(Ttf,g)] 



<\{Ttf,9)\ 



for all f,g £ ^^(Z^). Furthermore, it is easy to see that Ttf has an explicit representation: 



(5.24) 
where 



Ttf = f*{ m 



r(0) 



l(i7r)+iff); 



+ f*{\{H, 



(1) 



H 



{-!)■ 



(5.25) 



h[~'\x) 



H?\.) 



Hi'\x) 



1 



(2^)2 

1 

(2^ 
1 



Im 



Re 



-Im 



1 



J.{k-x-t'y{k)) ^^ 



[-7r,7r]2 l{k) 

^i[k-x-t"/{k)) ^^ 



[-7r,7r]- 



-f{k) e^C^-^-'TC^)) dk 



(see [12] for the finite volume analog). Combining (j5.23p and (|5.24p . we find that 
(5.26) \\[rt{Wif)),W{g)]\\<Y,\f{x)\\g{y)\ Yl \Ht\^ - y)\ ■ 

^,y mG{-l,0,l} 

We can now state the results. The first one describes pairs of Weyl operators W{f) and W{g), 
such that the norm of the commutator [Tt(W{f)), W{g)] decays polynomially. 

Theorem 5.1. Fix the parameters u > and Ai,A2 > 0. Denote by Tt the harmonic dynamics 
defined as above on W(£^(Z^)). Then there exists a number C3 > 0, such that 



(5.27) 



[Tt{W{f)),W{g)]\\< mill 



g3||/||i||g| 

|i|3/4 



holds for all f,g£f{Z'^). 

Next, fix a positive number 6 and denote X = supp(/) and Y = supp((/). Then there exist 
numbers Ci = Ci{5) > 0, i = 1,2, such that 



(5.28) 



[rt{W{f)),W{g)]\\< mm 



C, 



i\\9\\i 



|t|5/6 



holds for all f,ge i\Z^) such that X -Y eZ^\ BtsitVs), and 



(5.29) 



[Tt{W{f)),Wig)]\\<mm 



Ci 



i\\9\\i 



holds for all f,g€ 1^{1?) such that X — Y ^1?\ Bts{t{V2 U V3)). Here Br{S) represents an open 
neighborhood of radius r of a set S and X — Y is a difference set: 



X-Y = {x-y:x£X,y£Y}. 
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Proof. From (j5.26p . we have 

).W(q)]\\ < llflhllolli 



(5.30) \\[n{W{f)),W{g)]\\<\\fh\\gh max J^ |i/^)(x-y)|. 



me{-l,0,l} 

By Corollary 12. 4| for each m G {—1,0, 1}, there exists a constant Cm > 0, such that 

/^ 

(5.31) iHi'^Hx - y)\ < —^, for all \t\ > 1 and all x,y G Z^. 

This proves (lOTl) with C3 = C_i + Co + Ci. 

Next, we show (j5.28|) . Again applying ()5.26p . we have 

(5.32) l|[r*(VF(/)),W^(<7)]||<||/||iN|i max E l^^^(^-y)l- 

me{-l,0,l} 

By the assumption of (|5.28p . x — y^ Bts{tVs) for all x € X and y gY, therefore each integral in the 
right-hand side of the above inequality can be estimated by one of the following: (j2.1ip . (j2.12p . or 
the result of Theorem 12.71 Thus, for each m £ {—1, 0, 1}, there exists a constant C^ = C!^{6) > 0, 
such that 

C" 

(5.33) {H^'^Hx -y)\< --^, for ah |t| > 1 and ah x, y G Z^ such that x-y ^ Bts{tVi). 

This proves (lOTp with C2 = C2{5) = CLi{5) + C'q{5) + C[[5). 

The proof of (j5.29p is similar. D 

The following result is a direct consequence of Theorem 12.61 and Theorem 12.71 

Theorem 5.2. Fix the parameters a; > and Ajl,A2 > 0. Denote by Tt the harmonic dynamics 
defined as above on yV{f{I?)). Let X = supp(/) and Y = supp(5r) for f,g £ 1^{Z'^). Then if 
X — Y £ tVo with dist(X — Y, tV\) > 5, then there exist constants C = C{6) and fi = fj,{6) such that 

\\[rt{W{f)),Wig)]\\ < C||/||i||<7||ie-'^'ii^*(^-^'*^^). 

Moreover, for every fi > there exist constants < v^ < ^(1 + 2\/Ai + A2 sinh(///2)) and C^ < 
uj + 2i/Ai + A2 cosh(;u/2) such that 

(5.34) \\[rt{W{f)),W{g)]\\ C^ < e-M(dist(x,y)-.,|t|)^ 
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